Fractional Variational Calculus with Classical 



and Combined Caputo Derivatives* 



Tatiana Odzijewicz 
tatianaoOua . pt 



1 



Agnieszka B. Malinowska^ 
abmal inowskaOua . pt 



Delfim F. M. Torres^ 
delfim@ua.pt 

^Department of Mathematics 
University of Aveiro 
3810-193 Aveiro, Portugal 

^Faculty of Computer Science 
Bialystok University of Technology 
15-351 Bialystok, Poland 



We give a proper fractional extension of the classical calculus of varia- 
tions by considering variational functionals with a Lagrangian depending 
on a combined Caputo fractional derivative and the classical derivative. 
Euler-Lagrange equations to the basic and isoperimetric problems are 
proved, as well as transversality conditions. 

Keywords: fractional derivatives; fractional variational analysis; isoperi- 
metric problems; natural boundary conditions; Euler-Lagrange equations. 

2010 Mathematics Subject Classification: 49K05; 49K21; 26A33; 
34A08. 

1 Introduction 

Fractional calculus (FC) is a generalization of (integer) differential calculus, in 
the sense it deals with derivatives of real or complex order. FC was born on 
30th September 1695. On that day, L'Hopital wrote a letter to Leibniz, where 
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he asked about Leibniz's notation of nth order derivatives of a hnear function. 
L'Hopital wanted to know the result for the derivative of order n = 1/2. Leibniz 
rephed that "one day, useful consequences will be drawn" and, in fact, his vision 
became a reality. The study of non-integer order derivatives rapidly became a 
very attractive subject to mathematicians, and many different forms of frac- 
tional (i.e., non- integer) derivative operators were introduced: the Grunwald- 
Letnikow, Riemann-Liouville, Caputo [15l [171 US] , and the more recent notions 
of Cresson [TU], Jumarie [TB], or Klimek [T5] . 

The calculus of variations with fractional derivatives was born in 1996 with 
the work of Riewe, to better describe nonconservative systems in mechanics 
P71 [5S] . It is a subject of strong current research due to its many applications 
in science and engineering, including mechanics, chemistry, biology, economics, 
and control theory (see, e.g., the recent papers [2| [5l [6l [71 1^ [T3 l fT9 l [5T | [23 l [24] ) Fl 

Following [3S] , we consider here that the highest derivative in the Lagrangian 
is of integer order. The main advantage of our formulation, with respect to 
the "pure" fractional approach adopted in the literature, is that the classical 
results of variational calculus can now be obtained as particular cases. We 
recall that the only possibility of obtaining the classical derivative y' from a 
fractional derivative D°'y, a G (0, 1), is to take the limit when a tends to one. 
However, in general such a limit does not exist [55] . Differently from [5S] , where 
the fractional problems are considered in the sense of Riemann-Liouville, we 
consider here combined Caputo derivatives D"'^ . The operator "-^Dy^ extends 
the Caputo fractional derivatives, and was introduced for the first time in |22) as 
a useful tool in the description of some nonconservative models and more general 
classes of variational problems. More precisely, we investigate here problems of 
the calculus of variations with integrands depending on the independent variable 
X, an unknown vector-function y, its integer order derivative y', and a fractional 
derivative ^ D"'^y given as a convex combination of the left Caputo fractional 
derivative of order a and the right Caputo fractional derivative of order (3. 

The paper is organized as follows. Section [2] presents the basic definitions 
and facts needed in the sequel. Our results are then stated and proved in 
Section [31 We discuss the fundamental concepts of a variational calculus such 
as the Euler-Lagrange equations for the basic (Theorem [T^ and isoperimctric 
(Theorem [TS]) problems, as well as transversality conditions (Theorem [TS]) . We 
end with an illustrative example of the results of the paper (Section [4]) . 

2 Preliminaries on Fractional Calculus 

In this section we present some basic necessary facts on fractional calculus. 
For more on the subject and applications, we refer the reader to the books 

dSllITllSS]. 

^The literature on fractional variational calculus is vast, and we do not try to provide here 
a comprehensive review on the subject. We give only some representative references from 
2010 and 2011. Other references can be found therein. 
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Definition 1 (Riemann-Liouville fractional integrals). Let f G Li {[a,b]) and 

< a < 1. The left Riemann-Liouville Fractional Integral (RLFI) of order a 
of a function f is defined by 

arj{x) Hx-tr-'frn, 

r(a) Ja 

and the right RLFI by 

.I^f{x) f\t-xr-'f{t)dt, 

for all X G [a, b] . 

Definition 2 (Left and right Riemann-Liouville ft-actional derivatives). The 
left Riemann-Liouville Fractional Derivative (RLFD) of order a of a function 
f , denoted by aD"f, is defined by 

X S [a,b\. Similarly, the right RLFD of order a of a function f, denoted by 
xD'^f, is defined by 

X e [a, b] . 

Definition 3 (Caputo fractional derivatives). Let f e AC {[a. b]), where AC ([a, b]) 
represents the space of absolutely continuous functions on the interval [a, b] . The 
left Caputo Fractional Derivative ( CFD ) is defined by 

X G [a, b], and the right CFD by 

MiM - Jl- (-if) M - fit - .r'^fmt, 

X G [a, b], where a is the order of the derivative. 

Theorem 4 (Fractional integration by parts [TT]). Let p > I, q > I, and 

1 + ^ < 1 + a. If g d Lp {[a,b]) and f € Lq ([a, b\), then the following formula 
for integration by parts holds: 

b b 
J g{x)arj{x)dx = J f{x)xT^g{x)dx. 

a a 
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Definition 5 (The combined fractional derivative ^ D"'' 



]). Let a, (3 e (0, 1) 



and 7 G [0, 1]. The combined fractional derivative operator ^ D"'^ is given by 



Remark 6. The combined fractional derivative coincides with the right CFD in 
the case 7 = 0, i.e., ^D^''^f{x) =^ D'^f[x). For 7 = 1 one gets the left CFD: 



D^^^fix) =^ D^fix). 

Forf= [/i,...,/jv] : ^b] 
we put 



, iV e N, and /, e AC {[a, b]), i = 1, 



,N, 



^D!^'^f{x) := [^D!^'^f,{x),...,''D^^^Mx)] . 

In tlie discussion to follow, we also need the following formula for fractional 
integrations by parts f^ : 



gixfD^'^fix)dx = / f{x)Dt^n{x)dx 



'lf{xUl--g{x) - (1 - 7) fixUl-f^gi^x) 



(1) 



where D^'^^ (1 - 7),^^! + 7,D^ 



3 Main Results 

Consider the following functional: 

h 

J{y)= j L {x, y{x), y'{x),^ D^'^yix)) dx, (2) 

a 

where x G [a, b] is the independent variable; y(a;) G is a real vector variable; 
y'(a;) e with y' the first derivative of y; ^D'^''^y{x) € stands for the 
combined fractional derivative of y evaluated in x; and L £ C^ ([a, b] x M^^; R) . 
Let D denote the set of all functions y : [a, b] — > R^ such that y' and ^D'^'^y 
exist and are continuous on the interval [a,b]. We endow D with the norm 

llylli.oo := max l|y(x)||+ max ||y'(a;)|| + max ||'^L>^''^y(a;)|| , 

a<x<b a<x<o a<x<b 

where ||-|| is a norm in R^. Along the work, we denote by diK, i — 1, . . . , M 
(M € N), the partial derivative of a function K : R*^ R with respect to its 
ith argument. Let A S R''. For simplicity of notation we introduce the operators 
[■]"'^ and A i-}"'^ defined by 

[yp^ (x) := {x,y{^)y{^),'' Df'y{x)), 

X {yr/ (^) := C^, y{^),y'{^),'' D^^^ix), Ai, . . . , A.) . 
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3.1 The Euler— Lagrange equation 

We begin with the foUowing problem of the fractional calculus of variations. 
Problem 1. Find a function y G D for which the functional ([2]), i.e., 

b 

J{y) = J L[yr/{x)dx, (3) 

a 

subject to given boundary conditions 

y{a) = y-, y{b) = y\ (4) 
if ^ff £ IR^, achieves a minimum. 

Definition 7 (Admissible function). A function y E T) that satisfies all the 
constraints of a problem is said to be admissible to that problem. The set of 
admissible functions is denoted by V. 

Remark 8. For Problem [7] the constraints mentioned in Definition are the 
boundary conditions 

We now define what is meant by minimum of ^ on 2?. 

Definition 9 (Local minimizer). A function y £ V is said to be a local mini- 
mizer to J on V if there exists some 5 > such that 

J{y)-J (y) < 

for all functions y £ V with \\y ~ < (5. 

Similarly to the classical calculus of variations, a necessary optimality con- 
dition to Problem [1] is based on the concept of variation. 

Definition 10 (First variation). The first variation of J" at y £ D in the 
direction h £ T) is defined by 

SJ{y;h):^\im =—J{y + eh) 

provided the limit exists. 

Definition 11 (Admissible variation). An admissible variation at y £ T> for J 
is a direction h G D, h 7^ 0, such that 

• SJ^ {y; h) exists; and 

• y + eh £ T> for all sufficiently small e. 

Theorem 12 (see, e.g., |30) ) . Let be a functional defined on T>. Suppose 
that y is a local minimizer to J on T>. Then bj {y; h) = for each admissible 
variation h at y. 
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We now state and prove the Euler-Lagrange equations for Problem [TJ 

Theorem 13. If y = (j/i, . . . , y^v) a local minimizer to Problem]^ then y 
satisfies the system of N Euler-Lagrange equations 

d,L [yf/ {x) - ^di,+.L [y]-^;^ (x) + i^fl",92iv+.i [y]"/ (x) = 0, (5) 
i = 2, . . . ,N + 1, for all X e [a, b] . 

Proof. Suppose that y is a solution to Problem [1] and let h be an arbitrary 
admissible variation for this problem, i.e., 

h,{a) ^ h,ib) = 0, ^ = l,...,iV. 

According with Theorem 1121 a necessary condition for y to be a minimizer is 
given by 

d 



-g-^Jiy + eh)l^,=0, 



that is, 



'N+l 



N+1 



^''^ ^y^i''^ {x)hi^i{x) + ^ On+iL [y]'^''^ {x) — h,^i{x) 



i=2 
N+l 



J2^2N+^L[yr/{x) {^D^'Ph,^,{x)) 



i=2 



(6) 



dx = 0. 



Using the integration by parts formulas, for the classical and '-^D"'f^ derivatives, 
in the second and third term of the integrand, we obtain 



.N+l r 



/ E [^'^ [y]"'" (^) - [2/];'" (^) +^ Dty,^+.L [yr/ (x) 



hi_i{x)dx 



N+l 



^ hi^i{x)dN+iL [y]'^''^ (x) 



i=2 



c=6 



N+l 



Y,h,^^{x).Jl-^d2N+.L[yr/ (x) 



4=2 



-(1-7) 



N+l 



Y,h,.,(x)all-^d2N+.L[yr/ (x) 



x—b 



= 0. 



(7) 



Since hi{a) = hi{b) ^ 0, i — 1, . . . , N , we get 



N+l 

E 



d.L [yf/ (x) - -d^+.L [y]^^^ {x) +^ i^fl^S^^+.L [yr/ {x) 



hi-i{x)dx = 0. 



Equalities ([5]) follow from the application of the fundamental lemma of the 
calculus of variations (see, e.g., \2Qi)- □ 
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When the Lagrangian L does not depend on fractional derivatives, then 
Theorem fT3l reduces to the classical result (see, e.g., [30]). The fractional Euler- 
Lagrange equations via Caputo derivatives that one can find in the literature, 
are also obtained as corollaries of Theorem [T31 The next result is obtained 
choosing a Lagrangian that does not depend on the classical derivatives. 

Corollary 14 (Theorem 6 of [22'). Let y = (j/i, . . . ^hn) be a local minimizer 
to problem 

J{y)= L{x,y{x),^D^^Py{x)) dx^mm 

J a 

y(«) = y", y{b)^y\ 

y"> y** G M.^ . Then, y satisfies the system of N fractional Euler-Lagrange 
equations 

d,L[y]{x) + DtyN+^L[y]ix) = 0, (8) 
i = 2,...N + 1, for allx e [a,b]. 

If one considers 7=1 (cf. Remark [6]) and = 1 in Corollary [Ml then ([8]) 
reduces to the well known Caputo fractional Euler-Lagrange equation: if y is a 
local minimizer to problem 

Jiy)^ L {x, y(x), ^D'^y{x)) dx — > min 

J a 

y{a) = ya, y{b) = yb, 
then y satisfies the fractional Euler-Lagrange equation 

d2L {x,y{x),^D^y{x)) + .:.D^d:,L {x,y{x),^D^yix)) = (9) 
for all X € [a, 6] (see, e.g., [T^V 

3.2 Transversality conditions 

Let I e {!,..., N}. Assume now that in Problem [T] the boundary conditions ([4|) 
are substituted by 

y(a) = y", y.ib) = 2/^ i = l, . . . , ^ for z =^ and yi{b) is free (10) 

or 

y(a)=y^ y.{b) ^ yl i = 1, . . . , N ior i ^ I, and yiib) < yl (11) 

Theorem 15. If y ^ (2/1, • • • 1 2/iv) *s « solution to Problem [1\ with either (|10p 
or (jlip as boundary conditions instead of then y satisfies the system of 
Euler-Lagrange equations ([5]). Moreover, under the boundary conditions (jlOII 
the extra transversality condition 

On+i+iL [y]"'*^ (x) + -fxIb~°'d2N+i+iL [y]"^'^ (x) 

-{l-^)j',-f^d2N+i+iL[yr/{x)] =0 (12) 

J x—b 
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holds; under the boundary conditions the extra transversality condition 

On+i+iL [y]"''^ (x) + iJl'°'d2N+i+iL [y]'^'^ {x) 

-{l-^Ulrf'd^M+i+iLlyr/ (x) 

holds, with taking place if yi{b) < y^. 

Proof. The fact that the system of Euler-Lagrange equations ([5]) is satisfied is 
a simple consequence of the proof of Theorem [T3] (one can always restrict to the 
subclass of functions h e D for which hi{a) = hi{h) = 0, i = 1, . . . ,N). Let 
us assume that the boundary conditions are p^ . Condition (IT2|) follows from 
Suppose now that the boundary conditions are ([TT|) . Then, there are two 
cases to consider, (i) If yiih) < j/f, then there are admissible neighboring paths 
with terminal value both above and below yi{b), so that hi{b) can take either 
sign. Therefore, the transversality condition is P^ . (ii) Let yi{b) = yf. In this 
case neighboring paths with terminal value yi < yi(b) are considered. Choose 
hi such that hi{b) > 0. Then, e < and the transversality condition, which has 
its root in the first order condition ([7]), must be changed to an inequality. We 
obtain □ 

When the Lagrangian does not depend on fractional derivatives, then the left 
hand side of ([T^ and (IT^ reduce to the classical expression 9jv+i+i-^ {x, y{x),'y' {x)) 
(for instance, when iV = 1 and y(a) is fixed with y{b) free, then we get the well 
known boundary condition d^L (6, y{b), y'{b)) — 0). In the particular case when 
the Lagrangian does not depend on the classical derivatives, 7 = 1, = 1, 
and we have boundary conditions (fTO]) , then one obtains from Theorem [15] the 
following result of [1 . 

Corollary 16 (cf. Theorem 1 of [T|). If y is a local minimizer to problem 
Jiy)= / L{x,y{x),'^^D^y{x))dx-^mm 

J a 

y{a) = ya {y{b) is free), 
then y satisfies the fractional Euler-Lagrange equation (j9|) . Moreover, 
[,/i-"53i(x,y(a;),^Z?^y(x))]^^, = 0. 

3.3 The isoperimetric problem 

We now consider the following problem of the calculus of variations. 

Problem 2. Minimize functional ([S]) subject to given boundary conditions (|4]) 
and r isoperimetric constraints 

b 

g= (y) = J [y];'^ {x)dx = , J = 1, • • • , r, (14) 



< (13) 
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where € ([a, b] x R^^;R) and £ R for j ^ 1, . . . ,r. 

Problems of the type of Problem [51 where some integrals are to be given 
a fixed value while another one is to be made a maximmxi or a minimum, are 
called isoperimetric problems. Such variational problems have found a broad 
class of important applications throughout the centuries, with numerous useful 
implications in astronomy, geometry, algebra, analysis, and engineering. For 
references and recent advancements on the subject, we refer the reader to [21 
m [m [5D]. Here, in order to obtain necessary optimality conditions for the 
combined fractional isoperimetric problem (Problem [2]) , we make use of the 
following theorem. 

Theorem 17 (see, e.g.. Theorem 2 of [M] on p. 91). Let J,Q^,...,Q'' be 
Junctionals defined in a neighborhood of y and having continuous first variations 
in this neighborhood. Suppose that y is a local minimizer to the isoperimetric 
problem given by ([3]), ^ and (jl4p . Assume that there are functions h^, . . . , h' G 
D such that 

A = iflki), o-ki ■= ^Q^iV'i h''), has maximal rank r. (15) 
Then there exist constants Ai, . . . , G M such that the functional 

r 

satisfies 

5F{y- h) = (16) 

for all h e D. 

Theorem 18. Let assumptions of Theorem \17\ hold. If y is a local minimizer to 
Problem\^ then y satisfies the system of N fractional Euler-Lagrange equations 

9.Fa {yr/ (^) - ^dN+.F, {yr/ {x) + Dt"-y^2N+^F, {y}^^^ {x) = 0, (17) 

i = 2,.. .,N + 1, for all X e [a, b], where function F : [a, b] x M^^ x M'^ ^ M is 
defined by 

Fx {yr/ ix) L [y]^'^ (x) - ^ \,G^ [y];''' (x). 

Proof. Under assumptions of Theorem II 71 the equation (jl6p is fulfilled for every 
h S D. Consider a function h such that h(a) = h(&) = 0. Then, 

= STiy;h)^ {y + eh) 

'N+l N+1 , 

J2 ^-^ iy}''/ {x)h,-i{x) + ^N+^F^ {y};^^ {x)j-K_,{x) 

1=2 1=2 



N+1 

J2 d2N+^F^ {yr/ {xfD^'^h,.,{x) 

1=2 



dx. 
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Using the classical and the integration by parts formula ((J), and applying the 
fundamental lemma of the calculus of variations in a similar way as in the proof 
of Theorem [131 we obtain dTT]). □ 

Suppose now, that constraints ([14]) are characterized by inequalities 
6 

Q{y)^ J G^y];^^ ix)dx < ^„ j = l,...,r. 
a 

In this case we can set 

b b 

J (g^ [yf/ (^) - ^) dx + J i^,{x)f dx = 0, 

a a 

J = 1, . . . , r, where <j)j have the same continuity properties as y^. Therefore, we 
obtain the following problem: minimize the functional 

b 

J{y)^ j L {xMx),y'ix),^ D^^^yix),^^)) dx, 

a 

where (j){x) = [(/)i{x), . . . , <pr{x)], subject to r isoperimetric constraints 

b 



dx = 0, j = 1, 



and boundary conditions Q . Assuming that assumptions of Theorem [TS] are 
satisfied, we conclude that there exist constants Xj G M., j — 1, . . . ,r, for which 
the system of N equations 

d,F {x, yix),y'ix),^ D^'^y{x), X,, . . . , A„ <^(a;)) 

'^-^N+^F {x, yix), y'ix),^ D^'^yix), Ai, . . . , A„ ^{x)) (18) 



+ D^^'^^d2N+rF {x, D^'^y{x),Xi,. . . , A,, 0(x)) = 0, 



°'^d2N+iF {x,y{x),y'{x),^ 
i = 2, . . . , iV + 1, where = i + E Aj (g^ - ^ + (^f) and 

A,(^,(x)=0, j-l,...,r, (19) 

hold for all x G [a, Note that it is enough to assume that the regularity 
condition (jlSp holds for the constraints which are active at the local minimizer 
y. Indeed, suppose that I < r constraints, say , . . . ,Q'- for simplicity, are active 
at the local minimizer y, and there are functions h^, . . . , G D such that the 
matrix B ~ (bkj), bkj 5Q^ (y; h-*), k^j — 1, . . . , Z < r has maximal rank I. 
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Since the inequality constraints 5'+^, ■ ■ ■ are inactive, the condition ([TO]) is 
trivially satisfied by taking Aj+i = . . . = = 0. On the other hand, since the 
inequality constraints , ■ ■ ■ ,0'' are active and satisfy the regularity condition 
at y, the conclusion that there exist constants Xj & R, j — 1, . . . ,r, such 
that (dH]) holds, follow from Theorem[THl Moreover, (IT^ is trivially satisfied for 
j = l,...,l. 

4 An Illustrative Example 

Let a G (0,1),A^ = 1,7 = 1, and ^ G M. Consider the following fractional 
isoperimetric problem: 

J{y) = (y'ix) + ^DZy{x)fdx min 

JO 

G{y) = /' (y'i^) + ^D^yix)) dx = e (20) 

JO 

2/(0) = , 2/(1) = ^ i?i-c. (- (1 - t)'"") (dt. 

In this problem we make use of the Mittag-Lcfflcr function Ea{z). We recall 
that the Mittag-Leffler function is defined by 

oo 1^ 
fc=0 ^ ^ 

This function appears naturally in the solution of fractional differential equa- 
tions, as a generalization of the exponential function 9 . Indeed, while a linear 
second order ordinary differential equation with constant coefficients presents 
an exponential function as solution, in the fractional case the Mittag-LefBer 
functions emerge |17, . 

In our example pO)) the function F of Theorem [18] is given by 

F{x, y, y', '^DZy, A) = (?/ + ^D^yf ~\{y' + ^D^y) . 
One can easily check that y such that 

2/(2:) = ^ (- (a; - t)^-") Cdi (21) 

• is not an extremal for Q; 

• satisfies y' + ^D'^^y = ^ (see, e.g., [HI p. 328, Example 5.24]). 
Moreover, (PT|) satisfies the Euler-Lagrange equations ([T7)) for A — 2£_, i.e., 

(2 (2/' + ^D^y) ~ 20 + .D? (2 {y' + ^D^y) 2^ = 0. 
We conclude that (|2T|) is an extremal for the isoperimetric problem (|20|) . 



11 



Remark 19. When a — >■ 1 the isoperimetric constraint is redundant with the 
boundary conditions, and the fractional isoperimetric problem (j20[) simplifies to 
the classical variational problem 



1 

2, 



J{y) = 4 / {y'(x)Ydx — ^ min 

•^0 (22) 

y(0) = 0, 

Our fractional extremal (j2ip gives y{x) — ^x for i — 1, . . . ^ N , which is exactly 
the minimizer of (|22p . 



Remark 20. Choose ^ = 1. When a — >■ one gets from (PO]) the classical 
isoperimetric problem 



J{y) ^ (y (^) + y{^)) dx — > mm 
Jo 

G{y) = I y{x)dx = - (23) 
1 







y(0) = y.(l)-l 

Our extremal (|21|) js t/ien reduced to the classical extremal y{x) = 1 — of 
the isoperimetric problem (I23p . 



Remark 21. Lef a = i. Then ^ gives the following fractional isoperimetric 
problem: 



1 / „ 1 X 2 



J{y) = (y'{x)+^Dly{x)^ dx~^mm 

G{y) = / (2/'(:r) + ^Dly{x)) dx = ^ (24) 



2/(0) =0, y(l)-c(erfc(l)e+ 
where erfc is i/ie complementary error function defined by 

erfc(z) = —= / exp[—t^)dt. 
T/ie extremal (|2ip /or t/ie particular fractional isoperimetric problem (1241) is 



2/(:r) - e f c-Grfc(V^) + ^ - 1 
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